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Abstract 

We continue our studies of spherically symmetric self-similar solutions in the SU(2) 
sigma model coupled to gravity. Using mixed numerical and analytical methods we 
show existence of an unstable periodic solution lying at the boundary between the 
basins of two generic attractors. 



1 Introduction 

This is the third paper in a series aimed at understanding the structure of self-similar 
spherically symmetric wave maps coupled to gravity. In the first two papers ^ |2] we 
showed that for small values of the coupling constant there exists a countable family of 
solutions that are analytic below the Cauchy horizon of the central singularity. In this paper 
we wish to elaborate on the analysis of a periodic self-similar solution whose existence was 
only briefly mentioned in pp. We tried to make this paper self-contained mathematically 
but we refer the reader to for the discussion of the physical background of the problem 
and to [3] for more on the role of self-similar solutions in gravitational collapse. 



2 Setup 

We showed in that under the assumptions of spherical symmetry, equivariance and self- 
similarity the Einstein equations coupled to the SU(2) sigma field reduce to the following 
system of autonomous ordinary differential equations for three functions W{x), A(x), and 
F(x) 

W' = -l + a(l- W 2 )F' 2 , (1) 

A' = -2aAWF' 2 , (2) 

, . /s/ sin(2F) . , 

[AF'y = ^A-f (3) 

subject to the constraint 

I- A-2asm 2 F + aAF' 2 {W 2 -I) = 0. (4) 



Physically the functions W and A parametrize the metric, the function F parametrizes 
the SU(2) sigma field and a is the dimensionless coupling constant (a = corresponds 



to no gravity). We are interested in solutions of equations (l)-(4) starting at say x = 
with the following initial conditions (as explained in pP these conditions ensure regularity 
of solutions at the past light cone of the singularity) 

F(0) = |, F'(0) = b, W(0)=1, A(0) = l-2a, (5) 

where b is a free parameter (since the system has reflection symmetry F — ► — F we may 
take b > without loss of generality). The value A(Q) follows from the constraint (4). In 
what follows we shall refer to solutions of equations (l)-(4) satisfying the initial conditions 
JSJ as 6-orbits. We showed in [2] that for a < |, 6-orbits exist locally and are analytic in 
b and x. 

It follows immediately from (l)-(4) that a 6-orbit can be continued as long as \W\ < 1 
(since then A is bounded away from zero). However, if W hits ±1 at some x, then the 
solution becomes singular. Now we shall show that generic 6-orbits become singular in 
finite time. Throughout the paper we assume that a < \. 

We first show that if b is small, then W tends to -1 for some finite xa- To see this, let 
/ = (F — 7r/2)/&. Then, in the limit 6^0, equations (l)-(3) reduce to 

W = -1 + a(l - W 2 )b 2 f' 2 -> -1, (6) 

A' = -2aAW b 2 f' 2 — > 0, (7) 
sin(26/) 2f 

[Af) ~ b(w*-i) w^i' (8) 

with the initial conditions 

/(0)=0, f(0) = l, W(0) = 1, A(0)=l-2a. (9) 

The limiting equations (6) and (7) are solved by W = 1 — x and A = 1 — 2a. Substituting 
these solutions into (8) we get the equation 

(l-2oOr + ^ = 0, (10) 

whose solution (which can be obtained in closed form) behaves as fix) ~ (1 — 2a) + 
(2 — a;) In (2 — x) for x — > 2. Thus, the term (1 — W 2 )f' 2 in (6) stays bounded so, by 
uniform convergence on compact intervals, the solutions with sufficiently small b will tend 
to W = -1. 

Next, we show that if 6 is large then W tends to +1 and A — > for some finite ig. 
This time we define the variables 

e = fe 2 x, fc(0 = 6 a (l-W(s)) J = - f )• (11) 

Then, in the limit b — > oo, equations (l)-(3) reduce to (where now prime is d/d£) 

ti = 1- ah[2- ^) s' 2 ^ l-2ahs' 2 , (12) 



A' = -2aA ( 1 - A") s ' 2 -2aAs' 2 , (13) 



6 2 

6/i(2 - /i/6 2 ) 



(^'r = *i: {2s ! b L - (14) 
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with the initial conditions 



h(0) = 0, A(0) = l-2a, s(0) = 0, s'(0) = 1. (15) 

It follows from (14) and (15) that, in the limit b — ► oo, As' = 1 — 2a. Plugging this into 
equations (12) and (13), and using (15), we get the limiting solution 

= (1 - 2a)jl^laJ , = -L v /T3 4^(1 -V^R). (16) 

This solution becomes singular at £ = l/4a. By uniform convergence on compact intervals, 
we conclude that for solutions of equations (l)-(4) with large b (and nonzero a), the function 
W(x) attains a minimum and then tends to 1 at some x — > xb, while the function A{x) 
drops to zero at xb- 

To summarize, &-orbits tend in finite time to W = — 1 if b is small, and to W = +1 if 
b is large. In what follows, we shall refer to these two kinds of solutions as to type A and 
type B orbits, respectively. We show next that the sets of type A and type B orbits are 
open. 

Lemma 1. If W(x) > and A(x) < 1/2 — a for some x > then the orbit is of type B, 
i.e., there is a finite x such that lim W(x) = 1. Moreover, lim A(x) = 0. 

X — >Xo X—>Xo 

Proof: Substituting equation (4) in (1) we get 

„ 1 - 2a sin 2 F „ 1 - 2a 
W' = -2 + >- 2 +— — ■ 

Thus, if A(x) < 1/2 — a then W'(x) > so if W(x) > then W remains positive. But then 
by equation (2) A decreases which implies by (17) that W' remains positive (bounded away 
from zero in fact) and hence W must hit +1 in finite time. To prove the second part of the 
lemma, note that by equations (1) and (2) we have (using the abbreviation V = 1 — W 2 ) 

l ' V = f. (18) 

Assume that A(x ) > 0. Then (V/A)(xq) = and since (V/A)(x) > for x < x we get a 
contradiction. Hence A(xq) = 0. 

Corollary. Type B orbits are open. 

Proof: If the 6o-orbit is of type B then by lemma 1, A(x, bo) < 1/2 — a and W(x, b ) > 
for some x > 0. Hence for nearby b we also have A(x, b) < 1/2 — a and W(x, b) > and 
thus, again by lemma 1, the 6-orbit is of type B. 

Proposition 1. Type A orbits are open. 

Proof: First, note that if the orbit is of type A and W(x) > then A(x) > 1/2 - a 
(since otherwise the orbit would be of type B by lemma 1). But for W < 0, by equation 
(2) A' > 0, hence A{x) > A(x a ) > 1/2 — a for x > x where x is the point at which 
W(x ) = 0. Thus, A > 1 — 2a for type A orbits. Now, let 6 -° r bit be of type A and 
consider a nearby 6-orbit. By continuity, there is a point x\ such that W(x\,b) is close 
to —1, W(xi,b) < 0, and A(x\,b) is greater than, say, 1/2 — a. First we show that such 
orbits have W'(x, b) < for all x > x\. To see this, notice that from equations (l)-(3) 

W" = -2aF'A- x (aAWF'^{W 2 - 1) - WAF' + sin(2F)) , 
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hence at the first zero of W'(x, b) after x\ we have 



W" 



= -^UAW ±2VaV sin 2F). (19) 

w=o AV ' K ' 



The numerator is negative because A > 1/2 — a and W is close to —1 while the denominator 
is always positive, hence W" < which is a contradiction. Thus, W'(x, b) < and 
lim W(x,b) exists (if the orbit stays in the region). Now we show that W(x2,b) = — 1 

x — >xi 

for any such orbit for some xi > x\. To prove this assume that lim W{x) = W > — 1. 

x — >oo 

Integrating equation (1) we get 

pX r-X 

I W'dx = x x -x + a I VF' 2 dx, (20) 

J X\ J X\ 

which gives a contradiction as x — > oo because the last integral in i|20|) is finite by equation 
(2) (remember that A is bounded from below). 

We now know that both type A and type B orbits are open so there must be orbits that 
are not type A or type B, that is, orbits that stay in W 2 < 1 for all x. Call these type C 
orbits. We note that type C orbits are defined for all x > since A(x) > 1/2 — a for all x. 

Proposition 2. For any type C orbit we have 

liminf W(x) < and limsup VF(x) > 0. (21) 

Proof: Suppose that there is an x\ such that W(x) < —L < for x > x\ (this is equivalent 
to limsupW^x) < 0). Then from equation f !18|l 

V\' 2W -L 
— = < < —L 

A) a - a - 

for x > x\ and hence (V/A) < for some X2 > x\ which is a contradiction since (V/A) > 
for all x. 

Similarly, suppose that there is an x\ such that W(x) > L > for x > x\ (this is 
equivalent to liminf W(x) > 0). From lemma 1 we have that A(x) > 1/2 — a for x > xi. 
Thus 

'V\' 2W 



A) A ^ 2L 

for x > xi which is a contradiction since by lemma 1 (V/A) < 2/(1 — 2a) for all x. 

From Proposition 2 we see that type C orbits must oscillate at infinity about W — 
(unless lim W = 0) . 

Once we know that type C orbits exist we turn to their numerical construction. Nu- 
merics indicates that the structure of type C orbits is rather complicated for large a . In 
this paper we restrict our attention to small values a < 0.42, where the structure is simple. 
Namely for each given a there is a single critical value b* (a) such that 6-orbits tend to 
the attractor A (resp. B) if b < b* (resp. b > b*) and the 6*-orbit is of type C. In other 
words, the 6*-orbit is a separatrix lying between two generic attractors A and B. In the 
next section we give numerical and analytical arguments that the 6*-orbit is asymptotically 
periodic. 
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3 Numerical solution 



A straightforward way to determine the critical value b* is to take two values 6^4 and bg 
leading to attractors A and B respectively and then fine-tune to b* by bisection. This 
procedure yields a pair of b that are within a distance e from b* (where e is limited by 
machine precision). Such marginally critical 6-orbits exhibit a transient periodic behavior 
before eventually escaping toward W = ±1 (sec Figures 1 and 2). 
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Figure 1: The function W(x) for two marginally critical &-orbits for a = 0.38: the type A solution 
with b — b* — e (dotted line) and the type B solution with b = b* + e (dashed-dotted line), where 
e = 1CP 17 . Superimposed (solid line) is the periodic solution constructed by the straddle-orbit 
method. 
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Figure 2: The projection on the (A,W) plane of the same solution as in Fig. 1. The periodic 
solution is seen as the unstable limit cycle. 



5 



This suggests that the system has an unstable periodic solution and the &*-orbit belongs 
to its basin of attraction. In other words, the value b* corresponds to the intersection of 
the line of initial data (W = 1, F = 0, F' = b) with the 2-dimensional stable manifold 
of the periodic solution. In fact, if we take any two points Pa and Pg in the phase space 
which lead to attractors A and B, respectively, and perform bisection we obtain the same 
asymptotically periodic solution. This indicates that the stable manifold of the periodic 
solution is the boundary between the basins of attractors A and B. 

Since the periodic solution is unstable and numerically it is impossible to set initial 
conditions exactly on the stable manifold, we cannot obtain too many cycles of the periodic 
solution. Although in our case this is not a serious difficulty because the positive Lyapunov 
exponents are not large (see Figure 7 below), we would like to remark in passing that using 
so called straddle-orbit method due to Battelino et al. |S] one can pursue the unstable 
periodic orbit in principle forever. This procedure, which can be viewed as a series of 
bisections, goes as follows. At initial time we choose two points Pa{x — 0) and Pb(x = 0) 
which lead to different attractors A and B and perform bisection until the distance between 
the iterates Pa(0) and Pb(0) is less a prescribed S. Next we integrate the equations 
numerically starting from the current -Pa(O) and -Pb(O) until the distance between the 
trajectories exceeds <5. When this happens at some time x we stop the integration, assign 
the points Pa{x) and Pb(x) as current representatives and repeat the bisection. Iterating 
this procedure one can progressively construct a trajectory staying within a distance 5 from 
the codimension one stable manifold. The numerical solutions obtained by this method 
are shown in Figures 3, 4, and 5. 
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Figure 3: The profiles of the periodic solution for a = 0.38. 
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Figure 4: The phase portrait of the periodic solution for a = 0.38. 




Figure 5: The phase portraits of periodic solutions for different values of the coupling constant a 
ranging from 0.01 to 0.42. Asa^O the loop shrinks to zero and F' — ► oo. 

Using the fact that non-critical &-orbits become singular in finite time we can easily com- 
pute the positive Lyapunov exponent A of the periodic solution. To this end consider a 
marginally critical 6-orbit with b = b* — e. Such an orbit approaches the periodic solution, 
stays close to it for some time and eventually escapes along the unstable manifold to crash 
at a point xa where W{xa) = — 1- Therefore we can write 

X A % approach T Aperiodic "T" X escape? 

(22) 

where x approach, Xperiodic, and x escape denote the lengths of respective intervals of evolution 
(we say that the solution "escapes" if its distance from the periodic attractor exceeds 
a prescribed value). During the periodic interval the distance between the &-orbit with 
b = b* — e and the periodic solution grows at the rate proportional to eexp(Ax), hence 
Xperiodic ~ (— 1/A)lne. This implies that the number of cycles n during this interval 
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behaves as n ~ (— l/AT)lne, where T is the period of the periodic solution. The length 
of the escape interval does not depend on the number of cycles but only on the phase of a 
cycle at which the escape from the periodic solution takes place, hence x escape ~ /(lne), 
where / is a periodic function with period AT. Summarizing, we have 

xa ~ — -r lne + /(lne) + const. (23) 
A 

The numerical verification of this formula is shown in Figure 6. Using l|23|l we calculated 
the dependence of A on the coupling constant a - the result is shown in Figure 7. 
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Figure 6: For a = 0.2, the locus of the point of crash xa is plotted as the function of the 
logarithmic distance from the critical value In e. The fit to the formula (23) gives A = 2.029. The 
period of the wiggles, corresponding to the function /(lne), is equal to 2.887 in agreement with 
the predicted value AT (where T = 1.418 was calculated independently from equation (34)) . 
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Figure 7: The positive Lyapunov exponent A of the periodic solution as a function of the coupling 
constant a. 
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4 Perturbation series 



In order to construct periodic solutions we consider equations (l)-(4) with initial conditions 

F(0) = 0, F'(0)=c, W(0)=0, A(0) = (l + c 2 )-\ (24) 

where c is a free parameter and the value A(0) follows from the constraint (4). We claim 
that for sufficiently small a there is a unique c such that F(T) = tt, F'(T) = F'(0),W(T) = 
W(0),A(T) = A(0) for some T > 0. Since the system is invariant under the shift F — *■ 
F + 7r, we call such solution periodic. Now we shall construct the periodic solution in a 
perturbative way using the Poincare-Lindstedt method 

We define the new variable y — ^= where lo is the unknown in advance frequency. We 

remark that the rescaling of the independent variable by y/a is essential in order to have a 
well-defined limit for a — > while the rescaling by lo is introduced for convenience in order 
to have the fixed period 2w. In terms of y equations (l)-(3) transform to (J3 = y/a) 

luW' = /?[— 1 + lu 2 (1 — W 2 )F' 2 ], (25) 
A' = —2u)/3AWF' 2 , (26) 

. 2 (AF'Y = (3 2S ^1. (27) 

and the constraint (4) becomes 

1 - A - 2f3 2 sin 2 F + uj 2 AF' 2 (W 2 - 1) = 0. (28) 

We consider these equations on the interval < y < 2n with the boundary conditions 

F(0) = 0, F(2tt) = 7T, W(0) = 0, A(0) = Ao, (29) 

where the value of the constant Aq follows from the constraint (28). We seek solutions in 
the form of a power series in (3 



oo 

W(y,(3) =Y,P k W k {y), A(y,(3) = ^/3 fc A fc (y), F(y, 0) = ^ f3 k F k (y). (30) 

fc=0 fe=0 k=0 



The key idea of the Poincare-Lindstedt method is to expand the frequency in the power 
series 



oo 

W (/?) = ]T/?W (3i) 

k=0 



and to solve for the coefficients Uk by demanding that the solution contains no secular 
terms. Thus, we substitute (f^Uf) and (f5T]) into (25)- (28), group the terms according to 
powers of f3 and require that the coefficients of each power of /3 vanish separately. In the 
lowest order O(l) we get 

W (y) = 0, A (y)=(l + ^pj , F (y) = |, (32) 

where loq is yet undetermined. In the next order we get the equation ojqWI = (— 1 +w 2 /4), 
so to avoid a secular term we need to have loq — 2. Then, all 0(/3) terms are zero and in 
the order O(0 2 ) we get 

W 2 (y)=0, A 2 (y) = ~, F 2 (y) = i sm(y). (33) 
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Iterating this procedure with the help of Mathcmatica we calculated the perturbation series 
up to order 0((3 23 ). For example, up to order 0(/3 8 ) we have 

= 2 -Y + T-^ p8 + oin (34) 

W(y,P) = My) f + ^ P 5 + 25 = 5 7 6 (2y) + SiD(%) /? 7 + 0(n 05) 

= 1 , -2 + 4 cos(y) ^ | -4 - 8 cosQy) + 5 cos(2y) ^ 6 

2 2 8 16 

+ -60 + 204co S (y)-102cos(2y) + 52co S (3y) ^ | 10 
384 

F(y,/3) = | + ^M/? 2 + 16 sin(2,) /? 4 + 81 Sln(?/) = 21 ± Sin(3y) ffi + 

1656sin(y) + 900sin(2y)-616sin(3y) + 9sin(4y) o8 10 
+ ^ P +U\Jj ). {67) 

We recall that the "physical" frequency is equal to ^ so it diverges as (3 tends to zero 
(while the amplitude of oscillations goes to zero). In this sense the periodic solution is 
nonperturbative even though we constructed it by the perturbation technique. 

For small values of /3 the perturbation expansion converges quickly to the periodic 
solution constructed numerically (see Figure 8). As /3 grows the convergence becomes 
slower and we need to take many terms in the perturbation series to approximate well the 
numerical solution (see Figure 9). The fact that two independent ways of constructing the 
periodic solution agree, makes us feel confident that the periodic solution does in fact exist. 
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Figure 8: For a — 0.1 we plot the numerical periodic solution and superimpose the perturbation 
series (35). Even at this low order the agreement is very good. 
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5 Final remarks 

We showed above that for small values of the coupling constant a, the critical 6*(a)-orbit 
is asymptotically periodic as x —* oo. In the preceding papers [1,2] we showed that for a 
generic value of a, the &*(a)-orbit evolved backwards in x becomes singular as x — > — oo 
(which corresponds to the singularity at the center). However, there exist isolated values 
of a (called a n , n = 0, 1, ...) for which the 6*(a)-orbit is regular as x — > — oo. Combining 
this with the result obtained above, we conclude that for a finite set of isolated values a n 
(satisfying a n < 0.42) the Einstein-wave map equations admit self-similar solutions that 
are regular at the center and asymptotically periodic outside the past light cone. 
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